Search for Electron Neutrino Appearance in a 250 km Long-baseline Experiment 
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We present a search for electron neutrino appearance from accelerator produced muon neutrinos in 
tiie K2K long baseline neutrino experiment. One candidate event is found in the data corresponding 
to an exposure of 4.8 x 10^^ protons on target. The expected background in the absence of neutrino 
oscillations is estimated to be 2.4 ± 0.6 events and is dominated by mis-identification of events from 
neutral current n" production. We exclude the to Ue oscillations at 90% CL. for the effective 
mixing angle in 2-flavor approximation of sin^ 20^e(— | sin^ 2613) > 0.15 at Am^^ = 2.8 x 10~^eV^, 
the best fit value of the disappearance analysis in K2K. The most stringent limit of sin^ 26^e < 
0.09 is obtained at Am% = 6 x 10"^eV^ 

PACS numbers: PACS numbers: 14.60.Pq, 13.15.+g, 23.40.Bw, 95.55.Vj 



In 1998, the Super-Kamiokande (SK) collaboration re- 
ported evidence of neutrino oscillation based on atmo- 
spheric neutrino observations favoring large mixing be- 
tween and and a Ato^ near 2.2 x 10~'^eV^ 0. Sub- 
sequently, solar neutrino data from various experiments 
have indicated Ve disappearance as a result of neutrino os- 
cillations to other active neutrino flavors {vfj, or Uj) with 
large mixing and a Am^ near 5 x 10~^eV^ 0, y|- The 



KamLAND experiment also observes a deficit of reactor 
Ve consistent with the same parameter values as those 
in the solar neutrinos. Recently, the KEK to Kamioka 
long-baseline neutrino oscillation experiment (K2K) 
reported indications of — > Vx oscillation using an ac- 
celerator produced beam. The measurement of z^^ 
disappearance in K2K results in neutrino oscillation pa- 
rameters which are consistent with the values derived 
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from the atmospheric neutrino oscillations. 

Measurements of atmospheric and solar neutrinos sug- 
gest mixing between all neutrino flavors. The v^^ appear- 
ance is predicted at the same Am^ as the one measured 
using the atmospheric neutrinos (Am^j^^) in the frame- 
work of 3-flavor neutrino oscillations with certain param- 
eter values. The probability of i/g appearance in the 2- 
flavor approximation is expressed by 

P{v^, ^ V,) = sin^ 20^, sin^l .27 Am%L / E) , (1) 

where and Am^g are the effective mixing angle and 
mass squared difference between the mass eigenstates in- 
volved, respectively, for appearance; L is neutrino path 
length in kilometers; and E is the neutrino energy in 
GeV. In 3-flavor framework with "one mass scale dom- 
inance" 6], the effective mixing parameter is related to 
the mixing angle 613 by 

sin^ 29 ^ sin^ 923 sin^ 2^13 ~ ^ sin^ 2^13, (^23 - 7r/4). 

(2) 

In this letter, we report results from the first search for 
i>e appearance using beam in the K2K experiment 
sensitive to the Am^^^ region. 

In K2K, an almost pure (98%) i/^ beam with a mean 
energy of 1.3 GeV is produced with the KEK proton syn- 
chrotron (KEK-PS) 0. The fraction of is approxi- 
mately 1% and the remainder are v^. Twelve GeV pro- 
tons from the KEK-PS hit an aluminum target embedded 
inside a pulsed magnetic horn system which focuses pos- 
itively charged secondary particles, mainly pions, toward 
a far detector located 250 km far from KEK. The sec- 
ondary particles decay to produce a neutrino beam. The 
stability of the pulse-by-pulse beam direction is checked 
by monitoring muons from pion decay with a set of ion- 
ization chambers and silicon pad detectors following the 
beam dump. The measurements from these monitors 
show that the beam is directed to within 1 mrad of the 
far detector, SK which is a 50 kt Water Cherenkov 
detector located in Kamioka, Gifu Prefecture in Japan. 

To reject cosmic-ray and atmospheric neutrino back- 
ground, the global positioning system is employed at 
both the KEK and SK sites to synchronize between beam 
spills and events observed in SK 0. The neutrino flux 
at KEK is measured by a near detector complex con- 
sisting of a 1 kt water Gherenkov detector (IKT) and 
a fine-grained detector (FGD) system. The FGD con- 
sists of a scintillating fiber detector with segmented wa- 
ter targets (SciFi) 10'], a plastic scintillator hodoscope 
(PSH), a lead-glass calorimeter (LG), and a muon range 
detector (MRD) [ij. The flux at SK is estimated 
by extrapolating the measured flux at KEK using pre- 
dicted flux ratio between SK and KEK (far/near ratio). 
The far/near ratio is evaluated by a beam Monte Garlo 
(MG) simulation as a function of neutrino energy and 
is validated using secondary pion kinematic distributions 
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FIG. 1: Distributions of PID parameters for 32 single- 
ring events based on (a) Cherenkov ring pattern, and (b) 
Cherenkov opening angle. The distributions for the data 
(closed circles), oscillated MC with (sin^2S^^, Am^^j) = 
(1.0,2.8 X 10"^eV^) (solid histo grams) and expected Ue sig- 
nal with full mixing (dashed histograms) are shown. Shaded 
histograms are the NC component of i/^ MC. 

measured with a gas Gherenkov detector downstream 
of the horn system. 

Electron neutrino events in SK are selected assuming 
i^e charged current quasi-elastic (GG-QE) interactions in 
an oxygen nucleus, i.e. Ve + n ^ e + p. Since the proton 
momentum in the reaction is typically below Gherenkov 
threshold in water, only the electron is visible. Thus, a 
single electron-like Gherenkov ring is the signature of Ve 
appearance. The Ve contamination in the beam is esti- 
mated by the beam MG simulation which predicts the 
ratio of the number of i^e interactions to that of to 
be 0.9% in SK. This estimate is checked by a measure- 
ment of the i^e fraction at KEK using the FGD system. 
A more serious background comes from neutral current 
(NG) interactions where a single 7r° is produced and one 
gamma-ray from its decay is not reconstructed. The NG 
tt'^ production rate was measured using the IKT, con- 
straining the cross section. 

This analysis is based on the data taken between June 
1999 and July 2001, corresponding to 4.8 x 10^^ protons 
on target (POT). A total of 56 fully-contained events are 
obtained in the 22.5 kt fiducial volume of SK. Single- 
ring events are selected to enhance GG-QE interactions 
against the 7r° production background. The details of 
the selection criteria for fully-contained single-ring events 
are found in Ref. |^. Particle identification (PID) is ap- 
plied to reduce j/^-induced backgrounds. A Gherenkov 
ring produced by an electron is diffused by its electro- 
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TABLE I: Summary of the event reduction for v^, ap- 
pearance search at SK. The first column hsts each reduc- 
tion step, and the second gives the number of observed 
events after each selection. The numbers of expected back- 
ground from i/p and beam without neutrino oscillations 
are shown in the third and fourth column, respectively. 
The last column corresponds to the expected number of 
CC interaction events induced by Vg, oscillated from v^^ with 
(sin2 2e»^c, Am^J = (1.0,2.8 x IQ-^eV^). 





DATA 


Vfi w/o osc 


beam 


Ve from osc 


FCFV 


56 


80 


0.8 


28 


Single Ring 


32 


50 


0.5 


20 


PID (e-like) 


1 


2.9 


0.4 


18 


E„is >100MeV 


1 


2.6 


0.4 


18 


w/o decay-e 


1 


2.0 


0.4 
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magnetic shower and multiple scattering, while that pro- 
duced by a muon has a clear edge, and a low energy muon 
has a smaller opening angle than an electron. Both the 
Cherenkov ring pattern and opening angle are required 
to be consistent with an electron event. PID param- 
eters are calculated from a log-likelihood difference for 
the electron and the muon hypothesis. Distributions of 
these PID parameters are shown in Fig^ Negative val- 
ues of the parameters indicate an electron-like event. Dis- 
tributions of data are consistent with the oscillated v^^ 
MC with (sin226';,;,,Am2^) = (1.0,2.8 X lO-^eV^), the 
best-fit parameters of the disappearance analysis in 
K2K 5J, where 6*^^ is the effective mixing angle in 2- 
flavor approximation for j/^ — > oscillation. The visible 
energy is required to be larger than 100 MeV to reject 
low momentum charged pions and electrons from muon 
decays in which muons are below Cherenkov threshold. 
Finally, events which are followed by a decay electron sig- 
nal within a 30 /xsec time window are rejected. A small 
fraction of v^, interactions are also rejected by this cut 
when they are accompanied with decay electrons origi- 
nating from pions in inelastic interactions. The overall 
efficiency to select CC interactions from the oscillated v^. 
is 57 % for Am^ = 2.8 x IQ-^eV^. 

One event is selected as an electron candidate as sum- 
marized in Tabled While ring-counting algorithms eval- 
uate this event as a single-ring event, under a careful 
manual examination it reveals that the remaining PMT 
hits out of the reconstructed ring form an additional ring, 
and the invariant mass of these two rings is found to be 
consistent with the tt" mass. Thus, we conclude that the 
observed event is consistent with a 7r° background event 
from NC interaction. 

In our simulation to predict the number of background 
events, NC and CC pion production are modeled follow- 
ing Rein and Sehgal for the resonance region, and 
GRV94 Q with the correction of Bodek and Yang 



TABLE XL Systematic errors in the expected number of z/^ 
background in SK. 



Jun.1999 Nov.1999~Jul.2001 



NC cross-section 


+22% 
-27% 


+20% 
-25% 


Ring Counting 


+ 15% 
-13% 


+ 15% 
-13% 


Particle ID 


±11% 


±11% 


Energy Spectrum 


±14% 


±1% 


Far/Near ratio 


+ 15% 
-11% 


±6% 




±4% 


±4% 




±3% 


±3% 


POT normalization 


±0.9% 


±0.6% 


CC-nQE cross-section 


±1% 


±0.4% 


Total 


±36% 


+33% 
-31% 



for the deep inelastic scattering region. The axial vec- 
tor mass for QE and resonance pion production in the 
simulation are 1.1 GeV/c^ and 1.2 GeV/c^, respectively. 
The interaction models used in this analysis are the same 
as in the disappearance analysis in Ref. |5( except for 
the normalization of NC with respect to CC-QE cross- 
section. The reduction of the background events is sum- 
marized in Tabled 

The expected background from j/^ interactions in the 
case of no oscillation is estimated to be 2.0 events where 
the normalization is determined by extrapolating the ob- 
served number of events in the IKT 5]. It is estimated 
to be 1.9 events in the case of v^^ oscillation with 

(sin226', Am2) = (1.0,2.8 x IQ-^cV^). Since the back- 
ground is dominated by NC 7r° production (87%) and 
the oscillated Vr has the same NC interactions, it is in- 
sensitive to the disappearance oscillation parameters. 

The NC 7r° production cross-section in the MC sim- 
ulation is checked by a IKT measurement of tt*^ events. 
In the IKT, tt" events are selected by requiring two e- 
like rings whose invariant mass is between 85 MeV and 
215 MeV. Muon events are also collected by requiring a 
single //-like ring in the IKT as a reference to flux. 
The NC/CC-QE cross section ratio is calculated from 
the ratio of the number of tt'^ events to that of muon 
events. The 7r° sample in the IKT is dominated by NC 
interactions (87%), while the muon sample is dominated 
by CC interactions (97%). The measured NC/CC-QE 
ratio is consistent with the MC simulation; the ratios of 
data to MC is imt^'^^- To cover the allowed range, 
from 0.92 to 1.27, without changing the NC cross-section 
model in our MC, an uncertainty of 30% is assigned on 
the NC/CC-QE ratio. This error is used to estimate the 
systematic uncertainty in the background. The uncer- 
tainty in the cross-section ratio of CC interactions other 
than QE (CC-nQE) to CC-QE is estimated to be ±20% 
as in Ref. 

The uncertainty in the i^^-induced background is esti- 
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mated to be ±0.6 events for the 2.0 events. Contributions 
from various sources to the systematic error are summa- 
rized in Table UTI Since the horn current and target di- 
ameter were different in June 1999 from the other period, 
systematic errors are estimated separately for these two 
periods and properly weighted to obtain the total sys- 
tematic error. The uncertainty in the NC cross-section 
gives the largest contribution of 1l25%- Systematic er- 
rors from ring counting and PID are estimated by com- 
paring the shape of the MC and data likelihood distri- 
butions for cosmic-ray muons and atmospheric neutrino 
events. They are assigned to be and ±11%, re- 

spectively. Systematic errors from the neutrino energy 
spectrum (±1.0%) and far/near ratio (±6.0%) are esti- 
mated in the same manner as in Ref. [^. Systematic 
errors from the fiducial volume definition and detection 
threshold in the IKT {eiKr) and SK (esx ) are estimated 
to be ±4% and ±3%, respectively. 

The expected background from beam I'e interactions in 
SK is estimated to be 0.4 events, which is derived from 
the Ve/Vfj. flux ratio predicted by the beam MC simula- 
tion and the i^^ flux extrapolated from the IKT mea- 
surement. The systematic uncertainty in the number of 
beam Vf. events is estimated to be 0.11 events, which is 
dominated by the uncertainty in the Ve energy spectrum. 
The Ve/Vfj, ratio has been verified by a measurement of 
Ve events in the FGD . The Ve events in the FGD are 
selected by requiring 1) a vertex inside the SciFi fiducial 
volume, 2) an energy deposit in the PSH of greater than 
20 MeV, 2.5 times larger than expected from a muon, 3) 
an energy deposited in the LG of greater than 1 GeV, 
and 4) no corresponding hits in the MRD. During an 
exposure of 2.9 x 10^^ POT, 51 electron candidates are 
selected with an estimated background of 24 induced 
events. The v^jv^ interaction ratio is measured to be 
1.6 ± 0.4(siai.)lQ g(s?/s.)% which is in agreement with 
the beam MC prediction of 1.3%. 

The observation of one electron event in SK is consis- 
tent with the expected background of 2.4 events in the 
case of no oscillation. A constraint on neutrino oscilla- 
tions from Vft_ to is obtained by comparing the observed 
number of electron events with the expectation assuming 
oscillations. The expected number of electron events is 
calculated by 



osc 



+ nJ^^ + N, 



BG 



(3) 



where N^^'^ is the number of electron events induced by 
oscillated i^e, N^^^ is that induced by beam v^^ and N^^^ 
is that induced by both CC and NC interactions of 
and NC interactions of and Vr from oscillations. The 
— > Ve oscillation signal, N^^'~' , depends on the proba- 
bility of appearance expressed by Eq. ^ The number 
of beam Ve induced background, N^J-' — 0.4, is treated 
as a constant, since a contribution of Ve v^ oscilla- 
tion is negligible. The CC component of N^^'^ decreases 
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FIG. 2: The confidence interval for Ue. oscillations as 

a function of the effective Am^^ at 90% C.L. (solid line). 
Dashed line indicates 90% C.L. sensitivity of the experiment 
for the current statistics. The area to the right of each curve is 
excluded. Dotted hue shows the limit at 90% C.L. by CHOOZ 
assuming sin^ 26,ie = | sin^ 2Si3 . 



with v^^ disappearance observed in K2K and atmospheric 
neutrino experiments , depending on the survival proba- 
bility. 



P(i/^ ^ z/^) = 1 -sin^ 261^ 



,sin\l.27 ^ml^L/E), (4) 



where 9,,,, and Am^,,, 



are the effective mixing angle and 
mass squared difference for disappearance, respec- 
tively. In the present analysis, we assume 9^^ ~ j based 
on the nearly full mixing observed by atmospheric neu- 



trino experiments, and Am 



2 

fifi. 



"fie ' 



which is implied 



in the framework of 3- flavor neutrino mixing by the small 
mass difference found in solar neutrino experiments |l7j . 
Thus, Nexp reduces to a function of two parameters, 9f^e 
and Arn^g. 

A probability density function (PDF) for Nexp is con- 
structed from the Poisson distribution convoluted with 
the systematic uncertainty. Given the observation of 
one electron event, the systematic uncertainty has a very 
small effect on the derived confidence interval. The confi- 
dence interval of sin^ 20^e is calculated using the method 
suggested in Ref. In the calculation, the best-fit pa- 
rameters are searched for in the 2-dimensional parameter 
space with sin^ 29fj,e bounded in [0,1]. 

Figure |2 shows the limit on sin^ 20^e as a function of 
Am^g. The experimental limits on the neutrino mixing 
for the Ve oscillation hypothesis are given at 90% 

C.L. for a parameter region with Am^g > 6 x 10~* eV^. 
Neutrino oscillations from j/^ to Ve are excluded at 90% 



5 



C.L. for sin' 



> 0.15 at Aml^ 



2.8 X lO-^eV^. The 



most stringent limit of sin^ 20^e < 0.09 is set for Ato^^ = 
6 X 10"'^ eV^. The sensitivity of the experiment for the 
current statistics is also shown in Fig. |2| 

Assuming 3-flavor neutrino oscillations and CPT in- 
variance, our results can be compared to reactor ex- 
periments. CHOOZ has excluded sin^ 2^13 > 0.1 at 
Aml^ - 3 X lO-^eV^ 1^. This corresponds to a limit 
of sin^ 29f^e < 0.05 at Am^g ~ 3 x 10~^eV^ assuming 
^23 — 7i'/4, and consistent with the present analysis as 
shown in FigEl The limit on sin^ 26'i3 by CHOOZ is 
converted by assuming sin^ '26 — \ sin^ 26*13. 

The K2K experiment searched for —>■ i/^ oscilla- 
tions with accelerator-produced muon neutrinos travel- 
ing 250 km. This is the first experimental search for i/g 
appearance with sensitivity down to the Am^ suggested 
by atmospheric neutrino oscillations. A single electron 
candidate is found in SK. The observed event is consis- 
tent with the expected background event. The limit on 
the Ve appearance is obtained. At the best-fit parameter 
values of the K2K i/^ disappearance analysis, we set the 
90% confidence hmit of sin^ 29^,^ < 0.15. 
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